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Abstract 

We theoretically study discrete photonic lattices in more than three dimensions and point out 
that such systems can exist in continuous three-dimensional (3D) space. We study discrete diffrac- 
tion in the linear regime, and predict the existence of four-dimensional (4D) tesseract solitons in 
nonlinear 4D periodic photonic lattices. Finally, we propose a design towards a potential realization 
of such periodic 4D lattices in experiments. 



PACS numbers: 42.65.Tg, 42.70.Qs, 81.05.Xj 



Dimensionality is one of the key attributes of physical systems which determines their 
properties. This particularly holds for photonic lattices - linear or nonlinear photonic struc- 
tures conventionally considered in one- and two-dimensional geometry, where the flow of 
ight exhibits a plethora of intriguing phenomena which hold great potential for applications 
1, 2|. Paradigmatic phenomena occurring in nonlinear photonic lattices are discrete solitons 
1, 21. Their prediction in one-dimensional (ID) photonic lattices in 1988 (3| has awaited 10 

n n 

years for the first experimental observation [4J]. With the suggestion [5[ and the experimental 
realization of optically induced photonic lattices, both in one jfj] and two j3] spatial dimen- 
sions, novel excitations such as vortices js, 9] (that cannot occur in ID) were discovered and 
explored. One- and two- dimensional photonic lattices discussed here |l|, |2J are also referred 
to as waveguide arrays, because they are continuous along another dimension (second or 
third) along which the li ght propagates. Photonic crystals provide the opportunity to study 



versatile linear (e.g., see 



10]) and nonlinear phenomena [11] in three spatial dimensions. In 



those systems one usually applies (continuous) Maxwell equations, but discrete lattice mod- 
els may also be used when the so-called Tight-Binding Approximation (TBA) is applicable 
121 ] . However, to the best of our knowledge photonic lattices in more than three dimensions 
were not yet considered. Here we theoretically study discrete photonic lattices in more than 
three dimensions. We point out that such systems can exist in continuous 3D space, that 
is, their experimental realization is not hindered due to the properties of our space. The 
properties of discrete diffraction and four-dimensional (4D) tesseract solitons are presented 
in 4D linear and nonlinear periodic discrete lattices. Finally, we propose a design towards a 
potential realization of such periodic 4D lattices in experiments. 

The fact that a discrete lattice can have its dimension larger than the continuous space 



it is embedded in is known from complex networks [13|. For example, the Internet can be 
regarded as a network of dimension close to 4.5, and the network of airports close to 3, even 
though they are embedded on a 2D surface of Earth [l^]. Complex networks have been 
scarcely considered in optical systems. We point out at an interesting concept of complex 
networks of interacting fields called solitonets fl^ . Ilo] ]. where the interaction dynamics at 



each individual node in the system has infinite degrees of freedom Hi 1151- 

We begin with a paradigmatic model - the discrete nonlinear Schrodinger equation 



2 



(DNLS) - which describes dynamics of light in lattices of various dimensions: 

peN a 

where ip a is the complex amplitude describing the field, N a denotes the sites that are coupled 
to the site a, J is the coupling (hopping) parameter that we assume to be equal between 
all coupled sites, and 7 is the strength of the nonlinearity. Throughout the paper we use 



filly 

q6 



ly used to describe 
It is applicable 



the following normalization: XlcJ^al 2 = 1- This model was succesfu 
dynamics of light in ID and 2D photonic lattices (waveguide arrays) 
when the lattice wells are sufficiently deep, such that each well has a well defined resonance, 
and the coupling between different lattice sites is weak. Thus, one can think of this model 
as describing a system of weakly coupled high-Q resonators. 

In theory, any two pairs of resonators can be coupled thus yielding versatile structures of 
complex networks of resonators, which calls for a more rigorous definition of dimension. The 
dimensionality of any such network can be calculated by the following procedure: Let us 
choose one resonator and calculate the number of resonators [call it N(l)] that one can reach 



in I or fewer connections. The number N(l 
is somewhat altered from that used in 



scales as N(l) ~ l D when I — > 00. This procedure 
13j for usual complex networks due to the fact that 
it is the possibility of coupling rather than the Euclidean distance between the resonators 
that matters here. It is straightforward to verify that the dimensionality of a simple cubic 
lattice corresponds to half the number of nearest neighbors, but we emphasize that this is 
not a generally valid prescription. To see that note that the well known body centered cubic 
lattice (BCC) or face centered (FCC) lattices have dimension three as expected. 

Up to this point the theoretical model (CQ) was general in a sense that coupling between 
any two pairs of resonators was possible. From this point on we focus on 4D "simple cubic" 
lattices defined as follows: every resonator is labeled by four indices, a = (i,j, k, I), and it is 
coupled to eight resonators labeled by + k, I), (i — l,j,k,l), (i,j + l,k,l), (i,j — l,k,l), 
(i,j,k + l,l), (i, j, k — 1, 1), (i,j,k,l + 1), and (i,j,k,l — 1). We proceed with a discussion 
of light propagation phenomena encountered in these 4D " simple-cubic" lattices. The first 
question concerns diffraction of light, i.e., dynamics when nonlinearity is absent (7 = 0). The 
phenomenon of discrete diffraction has been addressed many times in ID and 2D systems 
(e.g., see 0, Q]), yielding a characteristic pattern of lobes spreading during propagation. In 
4D it is impossible to visualize such a pattern and therefore to characterize diffraction we 



utilize the concept of inverse participation ratio, I(t) = l/X)a-^(*)a) where 1(a) = |?/; a (t)| 2 
is the intensity of light. A typical question that we wish to address is: if we excite a 
single resonator, how does the excitation spread through the system? For a .D-dimensional 
"simple-cubic" lattice, propagation of the complex amplitude is given by 

v jd2 ... jl} - J] i°WJt), (2) 

a=ji—jn 

where the initially excited (at t — 0) site is «o = (0, • • • , 0) with amplitude 1, and J n is a 
Bessel function of order n. The evolution of inverse participation ratio (IPR) is asymptoti- 
cally then 



m J±X(2JtA ~(^) D , (3) 

\n=— oo / x ' 

171 ] to express the sum. In what follows, for 



where we have utilized a formula from Ref. 
simplicity we take J = 1. 

It is evident that the inverse participation ratio asymptotically increases as a power law 
with a slow logarithmic modulation, where the power-law exponent equals dimension. The 
logarithmic modulation would not have been present in continuous systems and its existence 
occurs because of the presence of the lobes. The effective number of sites excited by the 
light is smaller during diffraction in discrete systems in comparison to volume in continuous 
systems where the lobes are not present. Figure [T](a) illustrates diffraction IPR-dynamics 
for a finite amount of time in a "simple cubic" lattice in 1,2,3, and 4D. Superimposed 
on the results for an infinite lattice (red dashed lines), we also show dynamics for a finite 
size lattice (blue solid lines) with N = 2 24 sites and periodic boundary conditions. This is 
important because finite size effects will occur in practical realizations of 4D lattices. For 
example, our 4D lattice has the length of only 2 6 = 64 sites along one dimension and the 
effect is clearly visible for large enough times. 

Next we consider a 4D lattice of nonlinear resonators that we model with a DNLS Eq. 
(JTJ). A paradigmatic nonlinear phenomenon that occurs in nonlinear lattices are discrete or 
lattice solitons. Versatile types of 1D-3D discrete solitons were predicted and/or observed 
in optics including bright on-site solitons , staggered solitons p, ll8|, vortex solitons 



9[, and octopole solitons [19]. The simplest type of soliton that one can consider in 
4D lattices is a 4D on-site bright soliton that is centered on a single site. We have found 
this soliton by self-consistently numerically solving the stationary DNLS equation with the 




FIG. 1: (color online) Discrete diffraction in 4D lattices and tesseract solitons. (a) Evolution of 
the inverse participation ratio I(t) during discrete diffraction for linear lattices with D = 1,2,3,4 
(bottom up). Red dashed lines depict I(t) for an infinite lattice, whereas blue solid lines correspond 
to a finite size system with N = 2 24 sites and periodic boundary conditions. For D = 4 the 
boundary is reached at t ~ 15 and the IPR starts to oscillate, (b) A schematic illustration 
of a tesseract soliton. The phases at nearest neighbor nodes differ by tt, 7/J = 630. (c) Inverse 
participation ratio dynamics of a tesseract soliton with (black dotted line) and without (red dashed 
line) noise demonstrating stability (the two lines are on top of each other) . When the nonlinearity 
is turned off, the tesseract initial condition diffracts (solid blue). 



focusing nonlinearity 7 > (not shown). It occurs only above some threshold value of the 
nonlinearity (this also holds for 2D and 3D solitons 19|, |20( ) . The relevant parameter here is 
in fact the ratio 7/J, because DNLS can be scaled; this means that by reducing the coupling 
parameter J the effective nonlinearity 7/ J can be made stronger. On-site bright solitons 
are also found in 1D-3D photonic systems, and one can expect that other types of 1D-3D 
soliton excitations will also exist in 4D systems. 

Here we predict a novel type of soliton that occurs solely in 4D lattices: the tesseract 
solitons equally excite 16 sites of a 4D cube (i.e., tesseract), and therefore cannot exist in 



1D-3D systems. By using the methods outlined in Ref. [19j, we find that these solitons 
exist when the neighboring sites are tt out-of-phase as illustrated in Fig. [T^b). The stability 
of these solitons was checked numerically; Fig. [TJ(c) shows the inverse participation ratio 
dynamics of a stable tesseract soliton, and its discrete diffraction when propagated without 
nonlinearity present. Without the tt out-of-phase feature the intensity on the neighboring 
sites would not repel, and the excitation would collapse. 

Up to this point we have theoretically considered a model which can represent lattices 
in more than three dimensions, and analyzed some linear and nonlinear phenomena in such 
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FIG. 2: (color online) Illustration of the idea for creating photonic lattices in more than three 
dimensions via complex networks of optical resonators, (a) A schematic illustration of a (finite 
size) two-dimensional lattice of coupled resonators, (b) The coupling scheme which is topologically 
fully equivalent to the one in (a), despite the fact that resonators are located on a ID line. By the 
same token a 4D (discrete) lattice can be embedded in a continuous 3D space, see text for details. 



systems. However, one may say that these theoretical considerations are not more than 
academic curiosity because experiments are performed in 3D continuous space. We point out 
that discrete photonic lattices with dimensionality greater than three can exist in continuous 
3D space. To illustrate that fact compare a system of coupled resonators schematically 
illustrated in Figs. [2(a) and (b). (the coupled resonators are connected by lines). For 
the sake of the argument let us assume that the coupling parameters between all coupled 
resonators are equal, and zero otherwise. The configuration in Fig. [2(a) is evidently a 
discrete 2D lattice. However, the system sketched in Fig. [2(b) is fully equivalent to that 
of Fig. |2(a). Thus, a 2D network of resonators can be constructed by embedding these 
resonators in ID geometry (on a straight line), provided that connections between distant 
(in Euclidian sense) resonators can be made. By using this line of reasoning it is evident 
that the existence of discrete photonic lattices of dimension larger than three, describable 
with model (0Q), is not hindered by the dimensionality of our space. 

The limitations on the dimensionality and structure of experimentally realizable discrete 
lattices depend on our ability to construct coupled resonators that may be on distant loca- 
tions. This is not an easy task because usually high-Q resonators are coupled via evanescent 
coupling, which implies that they have to be close to each other because the light is tightly 
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FIG. 3: (color online) Proposal for the design of a 4D photonic lattice, (a) Illustration of a "simple 
cubic" 4D lattice. Waveguides enable coupling between distant resonators along the edges of 3D 
cubes. Evanescent coupling (hopping) along the diagonal enables construction of the 4th dimension 
for the discrete lattices, see encircled region enlarged in (d). For clarity only two sublattices are 
displayed in (a); the number of sublattices that can be added along the 4th dimension depends 
on the size of each individual cavity R, and the lattice size a\ of every 3D sublattice. See text for 
details. 



bound to the resonators and evanescent fields extend only a few wavelengths. However, in a 
recent work Sato et al. have demonstrated that two distant photonic high-Q cavities can be 



coupled by using an appropriate waveguide with mirrors at its ends 



16| . The realization of 



the cavities and coupler was made in a photonic crystal structure [16[. Importantly, many 
Rabi oscillations were observed in that system [l^ which provides a great promise for the 
future construction of 4D lattices. 

Two ingredients towards the realization of more than three-dimensional discrete photonic 
lattices are thus present: (i) the fact that more than three dimensional photonic lattices 
can in principle exist in our continuous 3D space, and (ii) distant high-Q cavities can be 
coupled as was demonstrated in Ref. 16[. A proposal of a completely specified photonic 
structure which would lead to an experimental realization of more than 3D photonic lattices 
is beyond the scope of the present work. However, here we propose a design which seems 
as a viable path towards the realization of 4D "simple-cubic" photonic lattices. It is based 
on a combination of waveguide coupling and evanescent coupling between resonators. The 
design is illustrated in Fig. |3j First consider a 3D simple-cubic 3D lattice of cavities coupled 
with waveguides. Let the unit vectors of this 3D simple-cubic lattice be ai, a2, and a3. The 
distance between two adjacent sites ai should be much larger than the size of an individual 
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cavity R. Consider now that we place two such simple cubic structures next to each other 
such that each site of the second lattice is displaced by the vector d = d(a.i + a2 + as)/ y/3 
from the first one, where d is slightly larger than 2R, so that evanescent coupling (i.e. 
hopping) from one sublattice to the other is possible [see Figs. [3]^a) and (b)]. By adjusting 
the distances d and a/ ^> R, the coupling constant between two sites of the same sublattice 
(coupled via waveguides) and the hopping parameter between two adjacent sublattices, can 
in principle be made approximately equal. By this procedure we have constructed a bi- 
cubic lattice. However, since a ; 3> d ~ 2R, we can add a number (roughly up to aiy/Z/d) 
of sublattices along the direction d thereby creating a finite size 4D lattice. If a\ is a bit 
more than 10 times larger than d, then a finite size 4D simple cubic lattice with 10 4 sites 
can be constructed. For resonators made of nonlinear materials this construction yields 4D 
nonlinear photonic lattices. 

Let us comment on the finite size effects and scaling associated with realizations of these 
"simple-cubic" 4D lattices. From the design scheme presented in Fig. [3] we see that along 
the fourth discrete dimension the photonic lattice would have to be finite; its size along this 
direction, i.e., the number of sublattices that one can add, depends on the ratio of the cube 
edges ai and the distance d between the resonators along the fourth discrete dimension. Thus, 
to build a larger 4D lattice one needs larger a;. In practice all lattices are finite (including 
1D-3D) and if they are sufficiently large the finite size effects would be unimportant. Our 
calculations show that 4D lattices with 10 sites (10 sites along one dimension) can already 
exhibit 4D behavior for some phenomena like solitons. 

In conclusion we have studied discrete photonic lattices in more than three dimensions, 
and pointed out that such systems can exist in continuous 3D space. We have studied 
discrete diffraction in the linear regime, and predicted the existence of 4D tesseract solitons 
in nonlinear 4D periodic photonic lattices. These novel structures would open the way for 
investigating new optical phenomena that one does not encounter in usual 1D-3D systems, 
but could also provide us with better understanding of dimensionalities beyond 3D which is 
of fundamental importance. We envision the study of versatile novel types of solitons and 
instabilities in these systems including vortex like structures, gap solitons, surface states 
and surface solitons, incoherent light dynamics, and studies of quantum optical phenomena. 
These structures could also yield schemes and opportunities for creating novel optical devices. 

Note added. When this work was in its final stages we became aware of a paper 21] 
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which proposes realization of a 4D quantum models by using ultracold atoms in optical 
traps, where the 4th dimension is encoded in the internal states of the atoms providing an 
extra degree of freedom. 

This work is supported by the Croatian Ministry of Science (Grant No. 119-0000000- 
1015). 
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